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INTRODUCTION 

This  study  considers  only  passive  vibration 
isolation  by  an  optimization  algorithm.  Most 
physical  structures  designed  for  dynamic  environ¬ 
ments  have  isolator  elements  to  attenuate  the  res¬ 
ponse.  Examples  of  engineering  problems  that 
could  benefit  from  this  work  are  ground  vehicle 
response  and  equipment  or  instrument  vibrational 
response.  Passive  damping  might  be  of  use  to  damp 
out  the  vibrations  of  large  space  structures  and 
these  algorithms  could  be  used  for  the  selection  of 
damping  parameters. 

The  constraints  considered  are  displacements  or 
accelerations.  The  frequency  constraint  has  not  been 
used  since 
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it  tends  to  destroy  the  linear  character  of  the  algorithm 
as  discussed  in  reference  (8) .  The  design  variables  are 
linear  changes  to  mass,  stiffness  or  damping  matrices. 
However,  only  numerical  examples  are  presented  for  stiff¬ 
ness  changes.  The  constraints  can  be  expressed  in  either 
the  time  or  frequency  domain  and  the  cumulative  constraint 
is  used  to  measure  the  amount  of  constraint  violation. 

The  objective  function  represents  a  design  variable  that 
restrains  displacements  or  accelerations  to  be  less  than 
a  maximum  value  at  a  single  response  point. 

It  is  shown  that  the  variation  of  displacement  or 
acceleration  constraints  are  shallow  in  reciprocal  design 
variables.  The  optimization  problem  formulated  in  this 
space  is  almost  linear.  The  weight  minimization  problem 
has  not  been  considered.  However,  the  form  of  the  dis¬ 
placement  or  acceleration  constraints  could  be  used  with 
weight  minimization  algorithms,  but  this  is  a  nonlinear 


optimization  problem. 

TJui  developed  algorithms  have  been  studied  for 
transient  response,  frequency  response  and  stationary 

random  using  the  direct  dynamic  solution.  The  algorithms  could  be  used 
with  a  reduced  basis  of  old  eigenvectors  as  well.  Multiple  response 

points  and  loading  conditions  may  be  used. 

TRANSIENT  RESPONSE 

The  minimization  of  displacements  or  accelerations  can 

be  formulated  as  a  MIN-MAX  optimization  problem  for  a  single 
•  # 

response  point  X^. 
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(1)  MIN  (MAX  |x\|  ) 

(2)  MX  +  CX  +  KX  «  P 

(3)  MAX  |X±  (t)  -  Xj(t)|-Xu 

(4)  K  «  K  +  lo.K. 

O  1  1 

(5)  M  *  M  +  Zo.M. 

O  11 

(6)  C  «  C  +  la. C. 

o  11 


(7)  u 


Only  the  direct  method  of  solution  has  been  considered 
in  this  study. 

Equation  (1)  minimizes  the  maximum  acceleration  in  the 
time  domain.  The  objective  function  could  be  displacements 
instead  and  the  present  algorithms  could  also  be  used. 
Equation  (2)  is  the  structural  dynamic  equations  in  matrix 
form  which  describe  the  displacement  response  X(t).  Equa¬ 
tion  (3)  is  the  so  called  relative  displacement  or  rattle- 
space  constraint.  The  present  algorithms  can  include  this 
type  of  constraint  in  the  analysis.  However,  no  specific 
numerical  examples  are  presented  using  the  rattlespace 
constraint.  Equations  (4),  (5),  and  (6)  show  the  linear 


changes  to  the  stiffness,  mass  or  viscous  damping  matrix 


with  the  design  variables  a^.  The  design  variables  could 
contain  differing  sets  in  equations  (4) ,  (5) ,  and  (6) . 
Equation  (7)  lists  the  constraint  limits  on  the  design 
variables  a^. 


FREQUENCY  RESPONSE 

Sometimes,  it  is  convenient  to  solve  vibration  problems 
in  the  driving  frequency  oj  domain.  This  is  true  for  pro¬ 
blems  which  have  experimentally  available  results  for 
transfer  functions.  Also,  for  stationary  random  analysis, 
the  frequency  domain  transfer  function  must  be  determined. 
Equation  (2)  is  transformed  to  the  steady  state  frequency 
domain  by, 

(8)  X  -  RE{Xoelut>  ,  P  »  RE{Poelwt) 
where  RE:  denotes  real  part  of 

i  -  v=i — 

ml  driving  frequency 

XQ:  amplitude  of  harmonic  response 

PQ:  amplitude  of  harmonic  loading 

For  the  harmonic  substitution,  equation  (2)  becomes, 

(9)  (-w2M  +  iuC  +  K ) X  *  P^ 

o  o 

The  amplitudes  XQ,  PQ  are  complex  numbers.  Equation  (9)  may 
be  solved  repeatedly  for  XQ  given  PQ  and  w  using  complex 
arithmetic.  It  is  more  convenient  to  use  the  real  displace¬ 
ment  components  in  the  analysis.  The  method  of  reference  (1) 


is  used  to  work  with  the  real  and  imaginary  components  of 


XQ  -  U-iV 
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The  optimization  becomes  in  the  frequency  domain  for  one 
response  point  Ui ,  . 

(11)  MIN  (MAX  (u)2/""u2  7^ 2)“ 
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(13)  MAX  |  Xoi(u>)  -  X  .  (u>)  |*  Xy 


(14)  K  -  K  +Ia.K. 
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Equation  (11)  is  the  amplitude  of  steady  state  accleration 
at  one  point  and  equation  (12)  are  the  structural  dynamic 
equations  to  be  solved. 
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STATIONARY  RANDOM 

A  Frequency  Response  solution  is  first  analyzed  to  de¬ 
termine  the  transfer  function  H(u)  which  is  either  the  dis¬ 
placement  or  acceleration  at  a  response  point  of  interest. 
The  spectral  density  of  the  output  is  given  in  terms  of  the 
spectral  density  of  the  input  for  a  single  input/output 
system  is  given  in  reference  (2) . 

S0(u)  «  | H (to)  |  2SI(w) 

The  same  reference  also  lists  techniques  for  analyzing 
multiple  input/output  systems.  The  mean  square  value  can  be 
calculated  for  any  frequency  interval , 


f  S  (u)  db). 
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Various  performance  measures  have  been  proposed  for 
random  analysis  such  as  using  either  the  spectral  density 
or  mean  square  value.  The  optimization  problem  for  sta¬ 
tionary  random  becomes, 

(18)  MIN  (MAX  S  ) 
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(20)  SQ  *  |H(w)  rsx 


(21)  K  ■  K  +  la.K. 
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The  maximum  displacement  or  acceleration  spectral  den¬ 
sity  is  the  objective  function  to  be  minimized  in  equation 
(18).  Only  the  single  input/output  case  is  used  to  calculate 
the  spectral  density  by  equation  (20)*  The  objective  function 
is  converted  to  a  set  of  equivalent  integral  constraints  and 
the  minimization  is  done  then  on  the  mean  square  response  in 
effect.  However,  the  algorithm  can  also  be  used  for  multiple 
input/output  with  minor  modifications. 

MIN-MAX  PROBLEM 

The  objective  function  (1) ,  (11)  or  (18)  can  be  converted 
to  a  simplier  algebraic  form.  Consider  equation  (1) , 

MIN (MAX | X ( t)  | ) 

This  minimization  is  equivalent  to  minimizing  an  addi¬ 
tional  design  variable  a  such  that 

MIN  a 

fxV  (t)  |  -  a  -  0  for  all  t. 

The  cumulative  constraint  has  been  used  in  the  optimal 
control  literature  (3)  to  convert  many  discrete  constraints 
in  the  time  domain  into  one  equivalent  integral  constraint 
which  measures  the  total  amount  of  constraint  violation.  In 
terms  of  the  cumulative  constraint,  the  MIN-MAX  part  of  the 


optimization  becomes. 


MIN  a 

(25)  /  <  IX^t)  |-  a  >  dt  -  0 

The  objective  functions  (11)  or  (18)  in  the  frequency 
domain  are  computed  in  the  same  manner  with  frequency  re¬ 
placing  time  in  the  integral. 

The  inner  problem  or  the  maximization  in  this  research 
was  done  by  function  evaluation.  This  is  efficient  for  the 
transient  problem,  but  the  frequency  response  problem  re¬ 
quires  a  decomposition  for  each  driving  frequency  in  equa¬ 
tion  (12) .  It  would  be  required  to  reduce  the  basis  of 
equation  (12)  by  using  the  real  normal  modes  for  efficient 
solution  in  locating  the  maximum.  Reference  (1)  recommends 
performing  a  one  dimensional  search  on  the  variable  u  to 
locate  the  maximum.  This  one  dimensional  search  would  re¬ 
quire  several  initial  starting  points  to  insure  convergence 
to  the  maximum  of  the  nonlinear  problem  in  w. 

ANALYTIC  DERIVATIVES 

For  statically  determinate  structures,  stresses  and  de 
flections  are  proportional  to  design  variables  that  are 
linear  changes  to  stiffness  such  as  areas  of  rods  in  truss 
members.  For  indeterminate  structures,  this  is  only  an 
approximation.  It  was  investigated  in  references  (4,5)  and 
found  that  high  quality  explicit  expressions  for  stresses 
and  deflections  could  be  generated  using  a  first  order 
Taylor  series  expansion  in  reciprocal  design  variables. 


That  is,  the  design  variable  space  for  stress  and  static  de¬ 
flection  is  shallow  in  reciprocal  design  variables.  The 
linearized  Taylor  series  expansions  represent  lines  that  are 
very  good  approximations  to  the  exact  constraints.  The  ex¬ 
pansion  of  a  response  quantity  <j>  is  done  in  the  reciprocal 
design  variables  0^  of  the  direct  design  variables  . 

*  *  *o  +  Z  ff^i 

®i  =  — 
ox 

The  direct  solution  of  the  dynamic  response  equations  in 
the  time  domain  uses  an  efficient  implicit  equation  solver 
such  as  Newmark  integration. 

This  would  be  the  most  general  capability  for 
solution  of  the  dynamic  equations.  The  Newmark  integration 
equations  presented  in  reference  (6)  are  listed  for  one  set 
of  integration  parameters,  6  =  *5  and  fc.  Given  the  response 
at  ,  the  response  at  t2  can  be  calculated  from  the  following: 

(26)  MX  +  CX.  +  KX.  =  P 

2  z2  Z2  Z2 


(27)  KK 


K  + 


M  + 


(28)  PP 

z2 


+  C  -rf  X.  +  X. 
At  t,  t. 


(29)  KK  •  X.  -  PP. 

Z2  Z2 


(30)  X.  =  4  (X  -  X.  )  -  4  X.  -  X. 

2  I?  *2  1  K  1 


(31)  X  -  X  +  At  X*.  +  At  *X 

r2  rl  2  '*1  2  z2 


The  displacements  at  the  next  time  step  are  calculate 
by  equation  (29) .  The  matrix  KK  is  only  factored  when  t! 
time  step  At  changes.  The  acceleration  and  displacement 
are  recovered  by  equations  (30)  and  (31) .  The  derivativt  - 
of  the  response  quantities  are  found  by  differentiating 
equations  (26)  through  (31)  as  was  done  in  reference  (7) . 
This  is  the  pseudo  loads  teachnique.  The  derivatives  with 
respect  to  the  reciprocal  variables  are, 


3X 

3PP. 

(32) 

KK  Z2 

3KK  v  . 

3Bi 

*  3*i  *2 
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(33) 

3KK  _  3K 
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2  3C 

Using  this  technique,  the  derivatives  of  displacement, 
velocity  and  acceleration  must  be  calculated  and  saved  for 
all  degrees  of  freedom  in  the  finite  element  model  at  two 
neighboring  points  in  time.  The  KK  matrix  in  (32)  was  de¬ 
composed  in  the  response  calculations  and  would  not  be 
factored  again  in  this  step. 

The  pseudo  loads  technique  was  applied  to  the  structural 
equations  (12)  in  the  frequency  domain.  The  required  deriva¬ 
tives  are. 


(37) 
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LINEARIZED  CONSTRAINTS 

The  acceleration  cumulative  constraint  has  the  following 
first  order  Taylor  series  expansion  in  the  reciprocal  variables 


(38) 


/T  I  dt 
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*  0 


•• 

=  X 


•  • 

3X. 


1  ST?""  -  B,  for  X  -  B 


where 
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♦  •  3X.  ,, 

1  "  -Xo  -  1  36 j  56i  -  B,  for  X  -  6 


1=0  For  X  otherwise 


This  constraint  is  numerically  integrated  by  a  modified 
trapezoid  law  which  finds  those  response  points  above  a  line 
for  which  the  constraints  are  violated.  The  algorithm  inter¬ 
polates  to  find  the  points  where  the  actual  constraint  is 
violated. 

The  steady  state  acceleration  amplitude  in  the  frequency 
domain  is , 

(39)  A  s=  a)2  A 2  +  V? 


The  first  order  Taylor  series  expansion  for  the  accelera¬ 
tion  amplitude  magnitude  is, 

(40)  /“  I  dto  =  0 


where  I  =  A  + 
o 


Cl) 


J  2  2 

UT  +  VT 
3  3 


1  <aef  +  36T)S6i  *  8 


The  spectral  density  acceleration  is  calculated  like 
equation  (40) .  The  square  root  spectral  density  of  the  out¬ 
put  S  is  minimized, 

S  =  w2|HU) 

This  is  the  form  of  equation  (39)  and  equation  (40) 
would  be  modified  by  multiplying  w2  by  (u)  .  A  formula 
could  be  derived  for  multiple  sources  with  cross  correlation 
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s 


correlation  in  a  similar  manner. 

SEQUENTIAL  LINEAR  PROGRAMMING 
The  problem  considered  in  this  study  of  minimizing 
a  linear  design  variable  subject  to  constraints  on  displace¬ 
ments  or  accelerations  in  the  time  or  frequency  domain  is  an 
almost  linear  problem  in  reciprocal  design  space.  It  is  only 
natural  to  use  sequential  linear  programming  as  the  optimiza¬ 
tion  algorithm.  A  primal-dual  linear  program  which  is  listed 
in  reference  (10)  was  used  as  the  optimizer.  Sequential 
linear  programming  is  described  in  reference  (11). 


NUMERICAL  APPLICATIONS 
Transient  Response: 

The  model  of  reference  (1)  shown  in  figure  1  was  sub¬ 
jected  to  the  displacements  inputs  f^(t)  and  f ^ ( t )  shown  in 
figure  2.  This  model  represents  a  vehicle  running  over  a 
bump.  The  transient  step  size  used  was  .1  sec  and  39  time 
intervals  were  calculated.  The  five  springs  were  used  as  de¬ 
sign  variables  with  limits  shown  on  figure  1.  The  objective 
function  was  a  design  variable  which  represented  the  maximum 
acceleration  at  point  1  in  the  model  over  the  3.9,  sec  time 
of  response.  The  acceleration  constraints  were  made  active 


when  the  acceleration  bound  was  99%  of  the  maximum.  Figure 
3  presents  the  decrease  in  acceleration  at  point  1  in  the 
model  versus  the  required  number  of  structural  analyses. 

The  linear  program  uses  design  variables  that  are  changes 
from  a  reference  and  the  change  can  be  positive  or  negative 
which  requires  two  variables  be  subtracted  to  keep  all 
variables  positive.  So  the  total  number  of  design  variables 
used  in  the  linear  program  was  eleven.  Initially,  the  reci¬ 
procal  variables  were  constrained  by  a  move  limit  to  lie 
within  -  25%  of  the  initial  values.  Convergence  was  obtained 
at  iteration  three.  The  spring  rates  found  at  the  optimum 
were, 

k1  =  51.2  lb/in 
k2  *  200. lb/in 
k^  =  200. lb/in 
k4  =  1600. lb/in 
k5  *  1000. lb/in 

2 

The  minimum  acceleration  obtained  was  228.8  in/sec  . 
Figure  4  presents  the  initial  response  versus  the  optimal 
one. 

Frequency  Responses 

The  model  shown  in  figure  1  was  subjected  to  equal  in- 
phase  displacement  inputs  f^t)  -  *2(t)  *  5  COSu)t  at  the 


tire.  This  would  represent  a  vehicle  on  a  shaker  table.  The 
five  springs  waj.e  used  as  design  variables  with  the  limits 
shown  of  figure  1.  The  acceleration  amplitudes  were  evaluated 
between  5  RAO/SEC  and  44  RAD/SEC  in  steps  of  1  RAO/SEC  as  the 
driving  frequencies.  The  objective  function  was  the  maximum 
steady  state  acceleration  amplitude  at  point  1  over  the  range 
of  driving  frequencies.  The  acceleration  constraints  were  made 
active  when  the  acceleration  was  99%  of  the  maximum.  Figure  6 
presents  the  decrease  in  acceleration  amplitude  at  point  1 
in  the  model  versus  the  required  number  of  structural  analyses. 
Initially,  the  reciprocal  variables  were  constrained  by  a  move 
limit  to  lie  within  -  25%  of  the  initial  values.  When  a  step 
was  not  minimizing,  the  percent  move  limit  on  the  reciprocal 
variables  was  decreased  by  50%  and  the  linear  program  was  re¬ 
solved  at  the  previous  design  point.  Convergence  was  achieved 
at  iteration  7  which  was  close  to  the  value  found  at  iteration 
5. 

The  spring  rates  at  the  optimum  were  as  follows: 

k1  =  52.9  LB/IN 

k2  *  231.2  LB/IN 

k3  =  215.1  LB/IN 

k4  =  1000. LB/IN 

k5  *  1311. LB/IN 


The  minimum  acceleration  amplitude  was  found  to  be 
2 

318.7  IN/SEC  .  Figure  7  compares  the  initial  acceleration 
amplitude  in  the  frequency  domain  versus  the  optimized  res¬ 
ponse. 

Stationary  Random  Response: 

The  model  shown  in  figure  1  was  subjected  to  a  random 

displacement  at  the  tire  patches  as  discussed  in  reference 

(13)  with  parameters  that  correspond  to  a  smooth  highway. 

A  Frequency  Response  solution  is  first  completed  with  a  unit 

i  ( wt  -  4 ) 

harmonic  displacement  with  phase  lag  e  at  the  rear 

tire  with  phase  angle  $  »  ^ .  The  spectral  density  of  the 
output  in  terms  of  the  spectral  density  of  the  input  and 
transfer  function  is, 

SQ(u)  »  |h(u>)  (  2  Sjfw). 

The  transfer  function  is  determined  by  using  the  accelera¬ 
tion  output  of  the  frequency  response  solution  due  to  the  unit 
harmonic  input.  The  spectral  density  acceleration  was  eva¬ 
luated  between  5  RAO/SEC  and  44  RAO/SEC  in  steps  of  1  RAD/SEC. 
The  objective  function  was  the  design  variable  representing  the 
maximum  acceleration  spectral  density. 

The  acceleration  constraint  was  made  active  when  it  was 
99%  of  the  maximum.  Figure  8  presents  the  decrease  in  the  ob¬ 
jective  function  versus  the  required  number  of  structural  ana¬ 
lyses.  The  reciprocal  spring  rates  were  constrained  by  -  25% 
of  the  current  value  as  move  limits.  At  the  detection  of  each 
infeasibility,  the  move  limit  was  reduced  by  50%  of  the  current 
percent. 


Convergence  was  achieved  at  iteration  3. 
The  spring  rates  at  the  optimum  were. 


Kx  *  51.2  LB/IN 
K2  *  200.  LB/IN 
K3  »  200.  LB/IN 
K4  «  1000.  LB/IN 
Kg  «  2000.  LB/IN 

2  2 

The  minimum  spectral  density  was  49.61  (IN/SEC  )  /Hz. 

The  initial  and  optimized  spectral  densities  are  presented 
in  figure  9  . 

REDUCED  BASIS 

The  dynamic  equations  are  usually  reduced  from  physical 
degrees  of  freedom  to  some  set  of  generalized  freedoms.  The 
following  transformation  is  used. 

X  =  YZ 

<YTMY)*Z  +  (YTCY)Z  +  (YTKY)Z  =  YTP 

Most  solutions  use  the  matrix  Y  as  the  collection  of  eigen¬ 
vectors.  When  the  equations  are  differentiated,  the  derivative 
of  the  eigenvector  must  be  calculated.  An  alternate  approach 
would  use  direct  Ri,tz  vectors  as  in  reference  (15).  When  small 
changes  are  made  to  a  structure  as  would  be  the  case  in  vibration 
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isolation  an  old  eigenspace  could  be  used  as  in  reference  (9). 

When  old  Ritz  vectors  are  used,  the  eigenvector  derivative 
is  not  calculated  and  the  previously  developed  algorithm  can  be 
applied  with  the  same  sequential  linearity.  Sparse  matrix  mul¬ 
tiplications  should  be  used  if  small  changes  are  made  AG  to  a 
matrix  G. 

ytg  y  +  ytagy 

o 

where  Gq  is  the  unchanged  part  of  the  matrix  and  is  calculated 
initially. 

LOCAL  MINIMA 

The  algorithms  presented  coverage  only  to  a  local  minima 
and  it  is  necessary  to  use  several  initial  designs  to  identify 
all  minima.  The  method  tends  to  converge  to  the  minima  which  is 
the  strongest  resonant  peak  that  is  closer  to  the  initial  design 
than  other  strong  peaks.  The  following  table  lists  the  initial 
design  and  optimized  with  the  minimum  transient  response  of 
figure  1. 

TABLE  I 

INITIAL  OPTIMUM  INITIAL  OPTIMUM  INITIAL  OPTIMUM 


200 

500 

500 

1800 

1800 


(LB/ IN)  Kx 

100 

K2 

300 

K3 

300 

K4 

1500 

Ks 

1500 

MIN  X  , 

( IN/SEu  ) 

22 

"a*  ***  “V/* 

238.4 


238.4 


MULTICRITERIA  OPTIMUM 


For  multiple  response  points  and  loading  conditions,  the 
techniques  of  multiobjective  optimization  is  required  and  multi¬ 
objective  programming  is  described  in  reference  (14).  The  simul¬ 
taneous  minimization  of  all  objectives  is  in  general  not  possible. 
Individual  optimization  is  done  on  one  objective  function  at  a 
time  with  the  remaining  objectives  treated  as  constraints  and 
bounds  determined  by  the  analyst.  A  multitude  of  solutions  are 
generated  depending  on  the  constraint  bounds  on  the  objectives. 

As  an  example,  the  minimization  of  response  for  the  model  was 
considered  using  all  of  the  three  previous  loading  conditions.  This 
problem  has  three  different  objective  functions  with  each  having 
different  units.  Each  objective  should  be  minimized  subject  to  con¬ 
straints  on  the  other  two.  To  illustrate  the  method,  the  transient 
response  was  minimized  subject  to  constraints  on  frequency  response 
and  stationary  random.  The  constraints  were  made  active  at  99%  of 
the  initial  design  or  the  minimum  of  the  maximum  response  of  any 
previous  iteration  during  the  sequential  linear  programming.  Conver¬ 
gence  was  achieved  at  iteration  6  which  was  very  close  to  the  results 
of  iteration  4.  The  results  can  be  presented  in  a  table. 


TABLE  II 

Iteration  4  Iteration  6 


K1 

102.8 

65.8 

K2 

200 

200 

K3 

200 

200 

K4 

1280 

1217 

K5 

1000 

1778 

Transient 

237.5 

235.7 

Frequency  Response 

377.2 

322.0 

Random 

49.4 

50.1 

An  infinite  number  of  solutions  can  be  generated  in  this  same 
manner  and  results  can  be  found  minimizing  the  other  objectives. 
Some  engineering  judgment  is  required  to  interpret  the  generated 
solutions . 

CONCLUSIONS 

An  almost  linear  optimization  problem  of  importance  in 
vibration  isolation  has  been  identified  and  algorithms  were 
developed  to  minimize  the  forced  vibrational  response  of  structural 
systems.  These  algorithms  should  replace  the  very  inefficient 
one  presented  in  reference  (9)  which  solves  a  series  of  practical 
problems.  The  linearity  depends  on  using  displacement  or  accelera¬ 
tion  as  the  only  constraints  in  the  time  or  frequency  domain.  The 
frequency  constraint  is  inherently  nonlinear  as  discussed  in  re¬ 
ference  (8)  and  it  has  not  been  considered  in  this  study. 

Only  the  direct  dynamic  solution  has  been  used,  but  a  reduced 
basis  of  old  eigenvectors  could  be  implemented  as  well.  Only 

local  convergence  has  been  shown  and  several  initial  design  points 
should  be  used  to  search  out  other  local  minima.  Multiple  response 
points  and  loading  conditions  may  be  used. 
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AMPLITUDE  XQ  *  5" 

VEHICLE  SPEED  S  =  450  IN/SEC 

WHEEL  BASE  L  *  120  IN 

d1  *  360",  d2  «  144* 
w.,  *  jr£  *  1.25TT 

w,  *  jrs  ■  3.125TT 

d2 

ti  *di/s*  *3  sec'  V  (d1+d2>s  =  1,12  SEC 

TIKE  LAP  FRONT  TO  REAR  tL 
*L  *  L/s  -  .2667  S€C 


FRONT  WHEEL  DISPLACEMENT 

f,  (*)  ■  X0(l  -  coswt) 

0£Utf 

f,  (*)  *  Xq(1  ♦  COSW  (t-tj)) 

REAR  WHEEL  DISPLACEMENT 

f2(t)  *  f i  < t-t2) 

FIGURE  2:  TRANSIENT  DISPLACEMENT  INPUT  FOR  FIGURE  1 


FIGURE  7 :  FREQUENCY  RESPONSE  OF  FIGURE  1 


ITERATION  NUMBER 

FIGURE  S:  STATIONARY  RANDOM  RESPONSE  OF  FIGURE  1 
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